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基于 分 解 集成 的 航空 货运 需求 区 间 预 测 研究 
李 智 "， 白 军 成 ? 


(1. 兰州 交通 大 学 交通 运输 学 院 ， 甘肃 730070; 2. 西安 电子 科技 大 学 经 济 管理 学 院 ,西安 710071) 


摘 要: 航空 货运 是 国家 重要 的 战略 资源 ， 在 国内 及 国际 间 的 贸易 中 扮演 着 不 可 或 缺 的 角色 。 对 航空 货运 需求 进行 
的 科学 预测 是 航空 公司 制定 基础 设施 规划 和 总 体 投资 决策 的 重要 依据 。 针 对 航空 货运 量 数据 的 不 确定 性 ， 从 实际 需 
求 出 发 ， 文 章 引入 Bootstrap 方法 进行 不 确定 性 估计 ， 提 出 一 种 基于 分 解 集成 的 区 间 预 测 方法 。 具 体 来 说 ， 首 先 用 局 
部 加 权 回 归 的 时 间 序 列 分 解 (STL) 方 法 将 货运 需求 数据 分 解 ， 其 次 由 支持 向 量 回归 (SVR) 和 季节 自 回 归 综 合 移 动 平均 
(SARIMA) 分 别 对 趋势 分 量 与 季节 分 量 做 预测 。 再 次 ， 创 新 性 地 将 白 噪 声 分 量 进行 提取 并 用 Bootstrap 方法 做 重 采 样 
处 理 。 最 后 ， 将 预测 结果 与 处 理 后 的 白 噪 声 进行 集成 重 构 ， 利 用 分 位 数 构造 区 间 进 行 不 确定 性 量化 。 对 中 国 两 大 枢 
纽 机 场 货运 数据 的 实验 结果 表明 ， 构 建 的 区 间 能 够 有 效 地 结合 预测 结果 量化 不 确定 性 ， 为 区 间 预 测 提供 了 一 种 新 的 
研究 思路 。 
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Research on interval prediction of air cargo demand based on decomposition integration 


Li Zhi!', Bai Juncheng? 
(1. School of Transportation, Lanzhou Jiaotong University Gansu 730070, China; 2. School! of Economics & Management, 
Xidian University, Xi’?an 710071, China) 


Abstract: Air cargo is an important strategic resource of country and plays an indispensable role in domestic and international 
trade. Scientific forecasting of air cargo demand is an important basis for airlines to make infrastructure planning and overall 
investment decisions. Aiming at the uncertainty of air cargo volume data, this paper introduced Bootstrap method for 
uncertainty estimation and proposed an interval prediction method based on decomposition integration from the practical 
needs. Specifically, this paper decomposed the historical data by Seasonal and Trend Decomposition using Loess(STL) 
method firstly, then forecasted the trend and seasonal components by Support Vector Regression (SVR) and Seasonal 
Autoregressive Integrated Moving Average (SARIMA) , respectively. Thirdly, this paper extracted and resampled the white 
noise component by Bootstrap method. Finally, the prediction results were integrated and reconstructed with the processed 
white noise to quantify uncertainty using quantile construction intervals. The experimental results of cargo data from two hub 
airports in China show that the constructed interval can effectively quantify the uncertainty in combination with the predicted 


results, which provides a novel research idea for probabilistic interval prediction. 
Key words: decomposition integration; interval prediction; STL; SVR; Bootstrap 


0 引言 对 外 客运 需求 量 。 但 数据 的 非 平稳 性 特征 使 得 传统 的 预测 模 
E 型 很 难 进行 精准 建 模 ， 平 稳 序列 建 模 等 前 提 假 设 则 进一步 限 
随 着 我 国 经 济 转向 高 质量 发 展 阶段 ， 电 子 商务 与 快递 物 。 制 了 其 应 用 名。 随 着 数据 挖掘 技术 的 兴起 ， 人 工 智 能 方法 以 
流 业 发 展 迅 速 ， 航 空 货运 需求 占 货运 市 场 的 比重 也 在 逐年 增 其 对 非 线性 数据 优越 的 预测 性 能 受到 了 大 量 的 关注 ， 被 广泛 
加 。 科 学 准确 地 预测 航空 货运 需求 ， 关 注 航空 货运 需求 未 来 ”应 用 于 各 预测 领域 。 邹 宗 民 等 人 包 为 预测 高 速 公 路 短 时 交通 
的 动态 趋势 ， 并 根据 其 动态 趋势 采取 相应 措施 ， 是 航空 运输 ，” 流 建 立 支 持 向 量 回 归 模 型 ， 结 果 相 较 反 向 传播 、 自 回归 移动 
业 持 续 健康 发 展 的 有 效 保障 ， 能 够 为 各 级 航空 运输 决策 部 门 ”综合 平均 模型 具有 更 低 的 误差 。 张 波 等 人 09 构 建 长 短期 记忆 
制定 发 展 战略 和 规划 提供 数据 支持 趾 。 然 而 ， 航 空 货运 量 数 。” ”神经 网 络 并 应 用 于 机 场 路 面 交 通 拥堵 预测 ， 分 析 认 为 该 模型 具 
据 自 身 的 高 噪声 、 不 确定 性 、 非 线性 和 非 平稳 性 特征 ， 使 得 。” 有 一 定 实用 价值 。 但 是 ， 不 管 是 传统 预测 模型 还 是 人 工 智能 模 
预测 极 有 挑战 性 。 因 此 ， 提 出 一 种 稳定 、 准 确 的 航空 货运 系 ”型 ， 单 一 模型 都 不 能 充分 反映 出 货运 需求 的 数据 特征 。 为 了 
统 预测 模型 是 一 项 有 价值 和 关键 的 任务 。 有 效 利用 各 种 模型 的 优点 ，Bates00 提 出 的 组 合 预测 思想 ， 受 
近年 来 ， 人们 对 货运 量 预 测 模型 进行 了 大 量 的 研究 ， 主 要 ”到 国内 外 越 来 越 多 研究 学 者 的 青睐 。 赵 建立 等 人 9 构建 卷 积 神 
包括 传统 预测 模型 、 人 工 智能 预测 模型 与 组 合 预测 模型 。 传 统 ”经 网 络 结合 残 差 网 络 的 模型 对 地 铁 站 进出 客流 量 进行 预测 ， 而 
预测 模型 ， 如 灰色 预测 法 P3、 回 归 分 析 法 内 和 系统 动力 学 方 ” ” 究 认为 此 模型 具有 比 传统 模型 更 好 的 预测 精度 。 赵 亚 斌 等 人 0 
法 回 ， 由 于 其 较 强 的 可 解释 性 与 对 线性 数据 预测 的 良好 效果 ， ”建立 改进 粒子 群 优 化 算法 与 相关 向 量 回 归 (MPSO-RVM) 组 合 
早期 在 货运 需求 预测 领域 得 到 普遍 运用 。 刘 笑 佟 等 人 史 结 合 不 ”模型 进行 短期 交通 流 预测 研究 ， 通 过 实验 证 实 了 模型 对 于 降 
确定 理论 建立 多 元 线性 回归 模型 预测 全 国 铁路 年 货运 量 。 张 。 噪 处 理 后 的 数据 有 较 好 的 预测 精度 。 除 了 对 原始 数据 直接 进 
亚 男 等 人 [构建 系统 动力 学 模型 , 预测 2020 一 2035 年 海口 市 。” 行 组 合 预测 以 外 ， 将 数据 进行 分 解 集成 组 合 预测 也 是 一 个 
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益 流行 的 思路 。 梁 小 珍 等 人 04 提 日 
运 分 解 集成 预测 模型 ， 选 择 双 优 
型 或 自 回归 移动 综合 平 坎 
良好 的 预测 精度 。Li 乞 
态 分 解 的 航空 货运 
oo 
优 于 基准 模型 的 精度 与 鲁 棒 性 。 
想 能 够 集成 各 个 模型 的 优势 ， 有 助 于 


的 问 


变化 > 
逐渐 受到 重视 


等 : 基于 分 解 集成 的 航空 货运 


于 奇异 谱 分 析 的 航空 客 


法 结合 支持 向 量 回 归 模 


elmer ei 
基于 变 分 模 态 分 解 与 经 验 模 
归 移 动 平 均 模型 与 布 
组 合 预测 ， 获 得 了 
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STL 时 间 序 列 分 解 
局 部 加 权 回 归 的 时 间 序 列 分 解 (Seasonal and Trend 


Decomposition using Loess, STL) 方 法 是 一 种 近邻 平滑 的 非 
参数 回归 模型 JJ。 在 迭代 过 程 中 ， 通 过 友 代 学 习 每 个 数据 点 


组 合 预测 ， 
英 型 预测 精度 不 


述 预 测 方法 在 


题 

上 上 |， 包含 的 信 
息 有 限 ， 难 以 表征 类 

信 

受 


性 。 为 了 获得 更 精确 的 数据 
和， 区 间 预 测 研究 
页 测 0617， 


负荷 预测 15， 


风速 进行 预测 。 


一 种 流行 的 方法 ， 如 


忆 


用 上 下 界 估 


计 法 (LUBE) 建 Y 


页 测 区 间 并 对 


法 ， 在 训 碳 多 的 三 


区 间 的 LUBE 方 
上 re 


产生 的 误差 进行 


的 鲁 棒 性 权 值 ， 在 微调 趋势 和 周期 性 方面 赋予 相对 较 小 的 权 
， 以 减少 离 群 值 的 影响 。STL 方法 中 的 局 部 加 权 回 归 过 程 
本 质 上 是 一 种 时 间 背 景 下 的 最 近邻 核 平滑 。 分 解 过 程 允许 
用 户 指定 趋势 和 季节 成 分 的 变化 量 ,也 允许 指定 周期 的 长 度 ， 
这 对 于 分 析 不 同时 间 分 辩 率 的 时 间 序列 是 有 用 的 。STL 方法 
有 三 个 重要 参数 : w -- 季 节 周 期 的 观测 次 数 ( 或 长 度 ); N, -- 季 


节 平 滑 中 局 部 加 权 回 归 平 滑 窗口 大 小 ， NN, -- 外 部 循环 鲁 棒 性 


的 迭代 次 数 。 其 他 参数 可 以 相应 地 确定 或 使 用 默认 设置 ， 如 
六 是 内 循环 的 迭代 次 数 ，X 是 趋势 平滑 中 局 部 加 权 回 归 的 平 
滑 窗口 大 小 ，WN 是 低 通 滤波 的 0 


列 分 解 为 趋势 分 量 、 季节 分 量 和 


We iit 


二 
以 


统计 研究 ， 提 出 高 斯 
预测 ，Serrano-Guerrero 等 
均值 与 标准 差 构 造 区 间 对 
用 不 同 分 布 函数 分 析 数 据 
测 区 间 的 构建 。 


测 扩 充 为 区 间 
芭 将 模式 识别 引入 预测 ， 
电力 负荷 进行 预测 ; 


ng 等 人 09] 利 
页 测 结果 实现 预 


现 有 ee 的 研究 取得 了 一 系列 


的 成 果 ， 但 


Nai 


还 是 停留 在 对 货运 


量 的 宏观 把 握 己 


更 需要 参考 短 


人 b) 缺 


其 他 领域 中 的 


少 关 于 货运 量 的 区 中 
以 制定 合理 的 运输 
求 波动 区 间 的 掌 


测 货运 需求 波动 区 间 
效 手 段 ， 对 货运 需 
风 险 应 对 策 各 。 


计划 ， 是 机 场 最 利 的 


利用 数据 


克 空 货运 需求 的 预测 研究 
期 变化 ed i 


全 国货 运 


， 决 策 者 


c) 鲜 有 文献 针对 数 


噪声 的 随 


机 性 影响 模型 处 理 
声 将 为 模型 预测 效果 带 来 极 大 提升， 


0 品 
车 于 以 上 问题 ， 本 文 从 


同 分 量 间 的 相互 影 
结果 可 以 用 加 法 模型 表示 : 


检验 P2。 构 造 的 统计 量 oo=ww+22Z2e-D， 其 中 X 是 样 
总 量 由 
为 4 的 卡 方 分 布 ，m 常 取 In(N) 。 给 定 显 著 性 水 平 < ， 则 拒 
域 是 0>x*ia4。 接受 原 假设 意味 着 原 序列 是 白 噪声 序列 , 否 
认为 序列 存在 相关 性 P33。LB 检验 就 是 为 了 验证 时 间 


密切 的 相关 关系 ， 历 史 数 据 对 未 来 的 发 展 有 一 定 影响 的 序列 ， 


， 从 而 使 预测 结果 更 加 准确 。STL 分 解 


所 


Y=T+S +R (1) 


其 中 ， 了 为 时 间 的 原始 时 间 序 列 ，T 、5, 、RR 分 别 表示 趋势 
分 量 、 季 节 分 量 和 剩余 分 量 。T +S 被 称 为 确定 性 或 可 预测 成 
。 去 掉 趋 势 和 季节 性 分 量 后， 剩余 分 量 中 的 随机 性 变 得 更 


明显。 


1.2 Ljung-Box 检验 (LB 检验 ) 


为 了 对 序列 的 随机 性 进行 推断 , Ljung 与 Box 提出 了 LB 


量 ， 是 样本 k 阶 沾 后 的 相关 系数 ， 该 统计 量 服 从 


序列 之 间 有 无 相关 关系 的 手法 。 只 有 那些 序列 值 之 间 具 有 


种 基于 


的 白 噪 声 分 量 


Bootstrap(BT) 方 法 的 航空 货运 测 方法 ， 将 数据 中 


,并 将 其 用 Bootstrap A 理 


结果 构建 预测 区 
Be ea 


闻 (CPIs), 以 此 降低 预测 难 
本 研究 的 创新 如 下 ， 


a) 提出 


后 结合 组 合 预测 
是 高 模型 精度 ， 


保 声 和 降 噪 后 数据 


才 值 得 挖掘 历史 数据 中 的 有 效 信息 ， 用 于 
1.3 ”SVR 模型 


页 测 未 来 的 发 展 中 |。 


沿 


支持 向 量 回归 (Support Vector Regression, SVR)P5 是 支持 


量 机 在 回归 领域 的 一 个 重要 应 用 ， 其 基本 思想 是 通过 非 线 


性 映射 函数 将 原始 数据 样本 集 映 射 到 高 维特 征 空间 已， 并 搜 


过 超生 奋进 行 线性 回归 。 在 接 下 来 的 步 又 中 ， 将 函数 以 一 种 


砚 测 区 间 的 思想 。 从 统计 分 析 


预测 二 果 构 造 
噪声 没有 分 析 价 


J 预 测 的 , 将 其 作 


度 而 言 ， 
含 任何 有 用 信息 ， 从 纯 随 
机 性 的 数据 是 不 
， 忆 此 ， 


预测 模型 在 数据 预 处 理 阶 段 将 


的 
白 
机 的 东西 中 找 不 到 
可 
本 


之 后 的 预测 


段 中 只 预测 名 


白 噪声 分 量 ， 
用 分 位 数 构造 


b) 提出 一 种 区 


网 辐 


度 。 对 于 
果 进 行 重 构 ， 再 


间 预 测 方法 。 该 方法 


生成 预测 区 间 ， 


保证 了 预测 的 可 靠 性 
及 任何 现 有 预测 方法 所 要 求 的 先 验 知 i 
设 ， 因 此 比 间 接 区 此 
c) 建立 STL、Bootst 
合 的 分 解 集成 模型 ， 
SARIMA- BT 模型 如 
解 后 的 趋势 分 量 


js 
二 


的 区 间 预 测 方 法 不 涉 
和 预测 误差 的 分 布 假 


rap SVR-SARIMA 组 合 预测 相 结 
区 间 的 构建 。 此 STL-SVR- 
寻 下 ， 根 据 STL 分 


础 上 结合 利 


预测 ， 利 用 各 模型 多 


用 合适 的 模型 进行 


并 在 此 基 


] Bootstrap ， 


为 建 区 间 。 


吕 


品 


验 风险 误差 。 C 越 大 ， 对 训练 过 程 中 误差 “ 越 大 的 样本 的 惩 
如 珊 越 大 。 
1.4 SARIMA 模型 


减少 复杂 性 的 方式 尽 可 能 地 简化 , 并 在 更 大 范围 内 进行 泛 化 。 
段 设 给 定 样 本 为 {o,y)=12… 相 ， 则 回归 形式 可 以 表示 为 


yC0 = CNCO+P 0O) 
PH ， 2 为 支持 向 量 ，2 为 偏 移 量 。 
引入 松弛 变量 5 和 可 以 较 好 地 解决 不 确定 性 ， 因 此 支 


持 向 量 回归 可 以 表示 为 


in 5lol + C2 -多 ) 


SLf 0) -YSe+ts, G) 
St.y— f(x)<eté, 
20207=12,...,l 


Pp，C 为 惩罚 因子 ，s 为 经 验 误 差 ，C 控制 SVR 模型 的 经 


委 
村 


作 自 回归 综合 移动 平均 (Seasonal Autoregressive 


Integrated Moving Average, SARIMA) 模 型 是 ARIMA 模型 的 


用 于 改进 ARIMA 模型 在 季节 时 间 序 列 建 模 中 的 性 能 。 


“ 展 , 
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ARIMA 模型 是 


著名 统计 学 家 Box 和 有 


线性 非 平稳 时 间 序列 模型 。ARIMA 的 
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局 数 据 库 (www.stats.gov.cn) 与 Wind 数据 库 


0 JenkinsP2g 推 广 的 一 类 
基本 原则 是 , 它 假定 过 


去 和 现在 的 情况 将 延续 到 未 来 。SARIMA(p, d, q)(P, D, Q)[S] 


模型 共有 七 个 参数 ， 分 别 是 非 季 节 自 回归 阶 数 ?， 季 节 自 下 


归 阶 数 P， 公 


式 。 正 则 和 季节 差分 算 


动 平均 部 分 的 ? 阶 和 


一 


差分 阶 数 4 ， 季 节 差 分 阶 数 D ， 非 季节 移动 
平均 阶 数 4 ， 季 节 移 动 平均 阶 数 2 ， 时 间 序 列 的 周期 9 。 模 
型 的 数学 表达 式 如 式 (4) 所 示 。 

a (B)A,(B)I- BYU-B)?Y =W+e (BEo(B) 
其 中 ， a,(8) 和 “%(3) 分 别 为 非 季 节 分 量 中 的 自 下 
4 阶 特征 多 项 式 。4(8') 和 Eo(8') 分 别 表 
示 季 节 分 量 中 自 回归 和 移动 平均 部 分 的 了 阶 和 40 阶 特征 多 项 
子 分 别 为 4-8) 和 4-B')。d 表示 不 同 


(4) 
日 部 分 和 移 


序列 的 次 数 ， 以 消除 趋势 对 序列 的 影响 ， 表示 不 同 序列 的 


次 数 ， 以 消除 季节 性 对 
为 固定 项 ，; 为 随机 误差 。 
1.5 ”Bootstrap 方法 


Bootstrap 方法 最 5 


讲 ， 它 是 一 种 基本 


| 学 家 Efron 提出 R71。 从 本 质 上 
和 蔡 换 重 采 样 的 非 参 数 数据 驱动 统计 


影响 。 为 1 时 刻 观 测 值 ， 久 


样 方法 ， 其 基本 思想 是 通过 对 原始 村 


随机 抽样 ， 构 造 


人 


这 些 自助 样本 对 


行 统计 推断 。 该 方法 只 


依赖 其 他 假设 ， 能 够 充分 控 


本 数据 进行 有 放 回 地 
总 体 分 布 进 


于 给 定 的 原始 观测 样本 数据 ， 不 
原始 观测 样本 数据 所 携带 的 总 体 


息 ， 适 应 性 比较 广泛 ， 在 小 样本 时 表现 出 的 优势 更 加 明显 。 


一 般 情 况 下 , 每 个 Bootstrap 数据 集 的 数据 点 的 数量 被 重 


新 取样 为 N 个 ， 使 其 
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Bootstrap 的 次 数 8 常 取 为 100、500、 
10000,， 并 且 已 知 当 8B 
于 计算 成 本 与 8 成 
时 间 ， 使 


型 。 模 型 可 分 为 四 个 模块 ， 即 数 
岂 与 区 间 构 造 模块 。 实 现 步 又 如 下 : 
1: 数据 分 解 模块 。 利 用 


曾 加 时 ， 可 以 准确 地 估计 
因此 有 必要 考虑 计算 机 性 能 和 


比例 增加 ， 
用 适当 的 3 值 。 


STL-SVR-SARIMA-BT 模型 


为 实现 对 于 航空 货运 需求 的 预测 , 考虑 到 其 非 线 性 特征 ， 
在 小 样本 数据 量 下 构建 基于 STL 


与 原始 数据 集 的 数量 相同 


。 所 使 用 的 


1000、2000、5000 和 


与 Bootstrap 的 


chinaXiv 


个 参数 。 


区 间 预 测 模 


昌 预 处 理 模 块 ， 预 测 模块 ， 


STL 将 原始 货运 数据 分 解 


为 趋势 分 量 、 季 节 分 量 和 剩余 分 量 。 由 于 不 同 机 场 数 据 间 的 


差异 性 ， 对 STL 模型 
步骤 2， 预 涡 
选择 的 


以 此 完善 模型 


的 参数 按照 数据 特 怕 
虽 各 分 量 的 特征 选择 预测 模型 ， 
性 。 选 取 能 预测 非 线 全 


分别 设置 。 


E 特 征 的 SVR 
E 效 果 显 著 的 SARIMA 模型 为 组 合 模 


到 的 趋势 分 量 与 季节 分 量 的 


型 输入 ， 分 


模型 与 处 理 周期 怕 
得 
别 


训练 SVR 模型 与 SARIMA 模型 。 


步骤 3: Bootstrap 模块 。 对 剩余 分 量 进 


预测 矩阵 。 


步骤 4: 区 间 构 造 模块 。 取 上 述 相似 预测 的 值 做 分 位 数 


上 述 模型 构成 所 用 方法 命名 模型 为 STL- 


处 理 构造 区 间 。 


明 其 为 白 噪声 序列 后 ， 对 


iu, 


Bootstrap 的 重 采 样机 制 对 剩余 分 量 随机 打 乱 ， 通 
机 影响 延续 模拟 未 来 的 随机 变化 。 处 理 后 的 剩余 


SVR- SARIMA-BT， 模 型 的 


3 ”实验 


3.1 数据 集 与 性 能 指标 


kk 体 框 架 如 图 


文章 选用 中 国 


国际 机 场 -- 北 京 机 场 (S1) 与 上 海 


机 场 (S2) 的 月 度 货 
2020 年 5 月 的 航空 货运 数据 


作为 下 


行 LB 检验 以 证 
其 进行 Bootstrap 处 理 。 利 用 
过 对 过 去 随 


| 练 集 部 分 作为 模 


别 加 到 预测 得 到 的 趋势 和 季节 分 量 的 总 和 中 ， 以 产生 相似 的 


1 所 示 。 


(www.wind.com.cn)。 文 章 在 预测 部 分 用 最 大 最 小 归 一 化 法 对 


数据 进行 处 理 


E， 归 一 化 公式 如 下 。 图 2 与 3 分 别 是 归 一 化 后 
几 场 月 度 货运 量 序列 折线 图 。 


y=(%, — Yrsin)/ Vrwx 一 yu) (5) 


STL 分 解 


“Tdens100g 


旭 1 预测 模型 框架 


Fig.1 Framework of predictive modelling 


二 述 性 过 


60 80 


时 间 


图 2 北京 机 场 
Fig.2 Monthly cargo volume series of Beijing airport 


1 
160 180 


FE。 从 1 


明度 可 看 日 


有 明显 的 不 对 称 性 。 此 外 ， 数 


计 ( 表 IJ) 可 知 , 月 度 货 运 量 数据 呈现 出 不 稳定 
上 ， 两 处 机 场 数据 的 1 
昌 的 峰 度 均 小 于 0， 说 


遍 度 绝对 值 均 大 


的 分 布 较 正 态 分 布 更 为 平坦 。 


为 了 定量 评 


有 60 8 10 1 1140 160 180 
时 间 


图 3 上 海 机 场 月 度 货运 量 序列 
Fig.3 Monthly cargo volume series of Shanghai airport 


十 模型 预测 性 能 ， 本 文选 取 三 个 统计 指标 : 


平均 绝对 误差 (M 
误差 (RMSE); 


履 盖 宽度 准 贝 
指标 ，MAE 月 


AFE)、 平均 绝对 百分比 误差 MAPE) 和 均 方 根 


个 区 间 评 价 指标 : 预测 区 间 覆 盖 率 (PICP)、 


预测 区 间 平 均 宽 度 (PINAW)、 平 均 覆 盖 率 误差 指数 (ACE) 和 
上 (CWC) 描 述 各 模型 性 能 。 有 具体 来 说 ， 对 于 统计 
于 直观 评价 模型 预测 精度 , MAPE 用 于 描述 预 
值 相 比 实际 值 的 偏离 程度 ，RMSE 对 预测 异常 值 敏感 ， 能 


究 对 象 ， 选 取 2006 年 1 月 至 
究 样 本 。 以 上 数据 来 自 统 


Ce 


LT 


恨 好 地 反映 出 预测 的 精密 度 。 对 于 区 间 评 价 指标 , PICP 用 


于 衡量 真实 值 落 在 预测 区 间 上 下 界 的 比率 , PINAW 描述 预测 
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区 间 的 狭窄 程度 ，ACE 反映 预测 区 间 与 区 间 名 义 置 信 水 平 判断 其 是 否 为 白 噪声 序列 。 针 对 两 个 白 噪 声 序列 的 LB 检验 
(PINC) 的 偏差 ,CWC 用 于 综合 评价 预测 区 间 的 覆盖 率 与 平均 。” 结果 如 图 6、7 所 示 。 
宽度 。 各 自 的 定义 如 表 2 所 示 。 由 下 图 可 知 ， 两 处 机 场 货运 量 数据 白 噪声 序列 LB 检验 1- 
对 于 分 量 预测 的 度量 ， y, ， 多 分 别 表示 t 时 刻 的 实际 值 12 阶 的 值 均 大 于 显著 水 平 0.05， 接 受 白 噪声 的 原 假设 ， 即 序列 
和 预测 值 ， 其 中 += 1, 2,…, N,N 为 样本 数 。 对 于 预测 性 能 指 ”” 均 不 存在 自 相 关 性 。 基 于 此 ， 可 以 说 不 将 剩余 分 量 纳 入 预测 过 
标 ，L 和 U, 是 预测 区 间 的 下 界 与 上 界 。R 是 目标 值 的 宽度 ， 程 ， 以 及 用 Bootstrap 打 乱 剩余 分 量 均 有 了 合理 的 前 提 。 
7 和 7 是 决定 PIs 性 能 的 两 个 超 参 数 。7 放大 了 PICP 的 差异 ， gs 
并 惩罚 无 效 的 PIs。7 为 控制 参数 ，4 为 区 间 置 信 水 平 。 值 得 站 一 一 一 ”~ 
注意 的 是 ， 当 PICP 高 于 预测 区 间 名 义 置信 水 平时 ，PIs 是 有 ”oo 
效 的 ， 因 此 建议 所 得 区 间 的 PICP 应 大 于 等 于 预先 确定 的 名 I 
义 置信 水 平 值 。 a 
将 数据 样本 分 别 以 常见 的 70%:30%,75%:25%,80%:20% TFT 
的 比例 划分 为 训练 集 与 测试 集 ， 由 预测 误差 最 小 原则 确定 最 CE 
终 划分 标准 , 预测 误差 结果 如 表 3, 可 知 将 数据 样本 以 75%:25% a 
划分 所 得 训练 误差 均 最 小 ， 因 此 将 数据 样本 据 此 比例 划分 得 站 we 
训练 集 样本 为 130 个 ， 测 试 集 样本 为 43 个 。 一 一 一 一 和 一 一 向 一 一 也 
表 1 机 场 货运 量 数据 基本 统计 描述 (万 吨 ) 灌 
Tab.1 Basic statistical description of airport cargo volume data 4 北京 机 场 货运 量 的 STL 分 解 结果 
(Ten thousand tons) Fig.4 Decomposition results of Beijing airport cargo volume by STL 
TT 机 场 ” 最 小 值 ”最 大 值 均值 标准 差 偏 度 峰 度 of 
© 北京 5.6000 16.1000 11.1785 2.7887 -0.1100 -1.1630 和 
C9 上 海 15.7200 39.3791 29.1513 5.0830 -0.3980 -0.5080 于 
© 表 2 评价 指标 定义 = a = Se 
[a Tab.2 Definition of evaluation index 台 册 
@ 指标 定义 表达 式 TYTTTFFTFTTTTFFFPTY 
ie MAE 平均 绝对 误差 MAE= > | -$I/N | 到 电 EE: 
名 MAPE 平均 绝对 百分比 误差 MAPE= > ly, 5) yl/N go A ea eA | 
AM RMSE 均 方 根 误差 RMSE = ND Gy, -$IN “00 
四 PiCP 。 预测 区 间 覆 六 率 Pen-,c /NG = sim 
= PINAW 预测 区 间 平 均 宽度 PINAW = 5 (UL)/NR 图 5 上 海 机 场 货运 量 的 STL 分解 结果 
2 ACE 平均 覆盖 率 误 差 指数 ACE = PICP_ PINC Fig.5 Decomposition results of Shanghai airport cargo volume by STL 
CWC ”覆盖 宽度 准则 CWC = PINAW (1+ ye ), y= | 下 
近 3.2 原始 数据 的 STL 分 解 10 
全 STL 分 解 阶段 的 目的 是 将 原始 数据 的 噪声 序列 提取 出 来 。 。 
O 噪声 的 存在 将 会 降低 数据 的 预测 精度 、 增 大 预测 的 难度 ， 而 呈 ? 
在 预测 之 前 将 噪声 剔 出 可 以 对 应 解决 上 述 问题 。 文 中 用 试 错 奔 6 
法 确定 模型 参数 ， 其 中 北京 机 场 货 运 量 数据 对 应 的 参数 ， 
(NMNN NAN ) 为 (12,1,5,13,21,17), 上 海 机 场 货运 量 数据 对 . 
应 的 参数 (和 Nj,NiNosNi,N,N; ) 为 (12,2,0,13,21,17)。 对 北京 、 上 海 2 
两 大 机 场 原 始 的 货运 量 数据 分 解 后 所 得 分 量 分 别 如 图 4 与 5 ge 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
所 示 。 可 以 看 出 ， 两 处 货运 量 数据 均 可 提取 到 明显 的 趋势 与 p 值 
季节 性 。 图 6 北京 机 场 白 噪声 序列 LB 检验 结果 
表 3 数据 集 划 分 标准 Fig.6 LB testresults of Beijing airport white noise sequence 
Tab.3 Dataset partitioning criteria 2 
数据 集 划 分 标准 MAE MAPE RMSE 
(训练 集 :测试 集 ) Sl S2 S1 S2 Sl S2 1 
70%:30% 0.0109 0.0048 0.0251 0.0081 0.0151 0.0067 , 
75%:25% 0.0038 0.0026 0.0057 0.00S2 0.0046 0.0044 注 7 
80%:20% 0.0190 0.0321 0.0304 0.0803 0.0238 0.0361 加 
将 分 解 得 到 的 趋势 分 量 ， 季 节 分 量 与 剩余 分 量 分 别 做 时 
间 序列 平稳 性 检验 (ADF)、 自 相关 函数 (ACF) 与 偏 自 相关 函数 3 
(PACF)， 分 析 各 分 量 的 平稳 性 与 自 相 关 性 。 由 于 篇 幅 所 限 ， ， 
下 面 仅 做 结果 说 明 o 综合 分 析 结 果 来 看 ， 趋势 分 量 与 剩余 分 0.00 0.20 0.40 0. 60 0.80 1.00 
量 均 为 平稳 性 序列 ， 趋 势 分 量 、 季 节 分 量 均 存 在 自 相关 性 。 bE 
剩余 分 量 的 ACF 与 PACF 可 初步 表明 其 不 存在 自 相关 性 。 图 7 上 海 机 场 白 噪声 序列 LB 检验 结果 


接 下 来 ， 对 分 解 后 得 到 的 剩余 分 量 做 LB 检验 ， 进 一 步 


Fig.7 LB testresults of Shanghai airport white noise sequence 
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3.3 SVR-SARIMA 预测 立 STL-SVR- SARIMA-GLS 模型 并 记 为 模型 S， 本 文章 所 提 


预测 阶段 只 将 原始 数据 的 非 噪声 序列 , 即 趋势 分 量 (Tb 与 出 的 模型 STL-SVR-SARIMA-BT 记 为 模型 6。 模 型 代号 与 简 
节 分 量 (S9 作 为 对 象 ， 由 于 剩余 分 量 已 被 证 明 是 白 噪声 , 将 和 所 示 。 为 直观 体现 模型 预测 效果 ， 对 预测 区 间 履 盖 


不 参与 预测 过 程 ， 而 是 另行 处 理 。SARIMA 模型 的 参数 设置 率 等 评价 指标 作出 详细 对 比 ， 对 比 结果 如 表 6、7 所 示 。 
为 : SARIMA(0,1,1)(1,2,0)[12]。 模 型 阶 数 的 设置 由 最 小 化 信 20 

息 量 准则 (Akaike Information Criterion,AIC)PS 与 贝 叶 斯 信息 1 

准则 (Bayesian Information Criterion，BIC)B29 经 过 多 次 模型 实 

验 所 确定 。SVR 模型 的 参数 设置 为 : 惩罚 因子 C 为 10, 经 验 

误差 6 为 0.01， 选用 多 项 式 核 函数 。SVR-SARIMA 模型 对 北 ol 3 R 

京 与 上 海 机 场 的 预测 结果 如 图 8、9 所 示 。 表 4 展示 了 单 模型 " 四 

与 组 合 模型 对 趋势 分 量 与 季节 分 量 的 预测 效果 对 比 ， 可 以 看 图 10 ”北京 机 场 货 运 量 区 间 预 测 结果 

出 S-S 模型 对 北京 机 场 与 上 海 机 场 货运 量 数据 预测 的 表现 均 Fig. 10 Forecast results of cargo volume interval of Beijing airport 
是 最 优 的 ， 其 误差 均 小 于 两 个 单 模 型 SVR 与 SARIMA 的 误 

差 ， 其 中 S-S 为 SVR-SARIMA 组 合 模型 的 简称 。 


en, 实 值 
- - LR-SARIMA 预 测 值 | 


图 11 上 海 机 场 货运 量 区 间 预 测 结果 


Fig. 11 Forecastresults of cargo volume interval of Shanghai airport 


0 10 20 30 40 
时 间 表 5 预测 模型 代号 与 简称 
图 8 ”北京 机 场 预测 结果 Tab.5 Code numbers and abbreviations for prediction models 
Fig.8 Prediction result of Beijing airport 模型 分 类 模型 代号 模型 简称 
40 | 可 和 2 
号 | 针对 预测 部 分 模型 2 STL-SARIMA-BT 
06 模型 3 STL-SVR-BT 
0 针对 噪声 处 理 部 分 模型 4 STL-SVR-SARIMA-MBB 
2 针对 区 间 构 造 部 分 模型 5 STL-SVR-SARIMA-GLS 
"I 二 计 所 提 模 型 模型 6 STL-SVR-SARIMA-BT 
生生 由 表 6、7 可 知 : a) 分 解 集成 组 合 预测 思想 在 数据 预测 中 有 
9 上海 机 场 预测 结 不 可 替代 的 作用 。 货 运 量 数据 中 包含 有 多 种 特征 成 分 ， 而 单一 
Fig.9 Prediction result of Shanghai airport 的 模型 仅 能 捕捉 部 分 特征 。 在 采用 了 分 解 策 略 的 模型 中 ， 本 文 
表 4 趋势 分 量 与 季节 分 量 预测 效果 采用 基于 分 解 集成 的 组 合 预 测 的 方法 与 其 他 单一 预测 模型 相 比 ， 
Tab.4 Prediction effect of trend component and seasonal component 有 更 高 的 PICP 值 ， 较 窜 的 PINAW 值 和 更 小 的 CWC 值 。 以 
分 量 MAE MAPE RMSE 北京 机 场 货运 量 数据 的 实验 结果 为 例 ， 在 置信 水 平 为 95% 条 件 
S1 S2 Sl S2 Sl S2 下 ， 所 提出 模型 的 CWC 值 均 比 模型 1，2，3 减 小 99% 以 上 ， 
THSI(SVR) 0.1339 0.0651 0.2451 0.1275 0.1790 0.0819 远 远 优 于 对 比 模型 。 由 此 可 知 , 单一 模型 预测 所 得 区 间 的 CWC 
TeHSt(SARIMA) 0.0803 0.0857 0.1046 0.1104 0.1181 0.1077 值 居 高 不 下 , 原因 既是 未 采用 “分 而 治之 ”策略 , 侧面 印证 了 组 合 


TerrSt(S-S) 0.0038 0.0026 0.0057 0.0052 0.0046 0.0044 预测 的 必要 性 。b) 用 Bootstrap 随机 打 乱 白 噪 声 剩 余 分 量 后 构造 
3.4 和 白 噪 声 的 Bootstrap 处 理 与 区 间 构 造 的 区 间 更 合理 。 以 上 海 机 场 货 运 量 数据 的 实验 结果 为 例 ， 在 

综合 考虑 处 理 效率 与 处 理 效果 , 本 研究 中 Bootstrap 的 次 言 水 平 为 80% 条 件 下 ， 模 型 6 的 CWC 值 比 模型 4 小 94.04%。 
数 8 取 值 1000。 将 剩余 分 量 ， 即 数据 的 白 噪声 部 分 随机 截取 ” ”分析 结论 产生 的 原因 , 模型 4 中 对 白 噪声 分 量 做 重 采样 所 用 
测试 集 长 度 ， 用 Bootstrap 方法 随机 打 乱 1000 次 ， 各 得 到 一 MBB 方法 是 根据 剩余 分 量 的 先后 顺序 将 其 划分 为 可 重印 
个 1000X43 的 矩阵 。 以 列 为 单位 对 数据 进行 由 小 至 大 的 排 ” 若干 个 块 ， 然 后 有 放 回 地 抽取 固定 长 度 的 块 中 的 所 有 数据 
序 , 之 后 各 列 取 选 定 的 分 位 数 (0.80,0.85,0.90,0.95) 得 出 区 间 的 ”成 新 的 样本 ， 适 合 处 理 存在 自 相 关 性 的 时 间 序 列 。 考 虑 到 


上 下 界 值 , 从 而 得 到 上 下 界 序列 ,图 10、11 分 别 显示 了 北京 ， 噪声 序列 不 存在 自 相关 , 用 Bootstrap 打 乱 白 噪声 序列 描述 
上 海 国际 机 场 货运 量 区 间 预 测 PINC 为 95% 的 预测 结果 。 据 中 的 不 确定 性 ， 以 此 模拟 未 来 序列 的 随机 成 分 更 为 合理 。 
3.5 ”对比 模 型 性 能 分 析 在 与 最 新 区 间 构 造 方法 的 比较 中 , 本 文 方法 的 PICP 更 接近 对 
为 了 进一步 验证 所 提 模 型 的 有 效 性 ， 本 文 建立 了 一 系列 ”的 PINC, 有 较 窗 的 PINAW 与 更 小 的 CWC。 一 般 情况 下 , 较 


对 比 模型 ， 其 按 类 型 可 分 为 针对 预测 部 分 、 针 对 噪声 处 理 部 。 的 PINAW 可 以 获得 理想 的 PICP， 但 本 文 期 望 PINAW 尽 可 
分 与 针对 区 间 构 造 部 分 的 模型 。 各 对 比 模型 的 构造 介绍 如 下 : 小 ， 而 CWC 可 达到 PINAW 和 PICP 间 的 权衡 。 由 评价 指标 
首先 ， 针 对 趋势 分 量 与 季节 分 量 的 预测 部 分 ， 设 计 选用 单一 ”可 看 出 所 提 模 型 6 的 CWC 值 在 各 置信 水 平 内 均 优 于 或 次 优 
的 线性 回归 (LR),SARIMA,SVR 模型 预测 趋势 与 季节 分 量 ， 他 模型 。 也 就 是 说 ， 比 起 模型 5 基于 点 预测 误差 将 确定 性 
对 应 模型 1， 2，3。 其 次 ,针对 噪声 处 理 方法 ,将 所 提出 框架 ”” 测 结果 扩充 为 区 间 ， 本 模型 构造 区 间 的 方法 可 得 到 包含 更 多 观 
中 的 Bootstrap 方法 用 Moving Block Bootstrap(MBB) 代 蔡 ， 测 值 且 宽度 更 小 的 预测 区 间 。d) 综 合 来 看 ,文章 提出 模型 的 PICP 
建立 STL-SVR- SARIMA-MBB 模型 并 记 为 模型 4。 最 后 , 针 ” 均 大 于 对 应 的 PINC， 即 预测 区 间 的 概率 是 满足 的 。 此 外 ， 模 型 
对 整个 区 间 构 造 部 分 ， 将 Wang 等 人 07 所 提出 的 以 预测 误差 ”的 PINAW 也 相对 较 小 ， CWC 值 在 各 置信 水 平 中 都 是 最 小 或 次 
分 布 结合 点 预测 结果 的 区 间 构 造 方法 应 用 于 本 研究 对 象 ， 建 小 , 这 表明 提出 的 STL-SVR-SARIMA-BT 模型 在 构造 短期 货运 
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需求 预测 区 间 方 面 比 其 他 方法 更 好 。 STL- SVR-SARIMA-BT 分 解 集成 组 合 模型 , 对 机 场 货运 需求 
表 6 区 间 预 测 结果 评价 (北京 机 场 ) 进行 预测 ， 同 时 引入 多 种 预测 模型 作为 基准 模型 进行 对 比分 
Tab.6 Evaluation of interval prediction results (Beijing airport) 析 ， 提 供 了 一 种 构造 区 间 的 方法 。 模 型 的 核心 思想 是 将 数据 
PINC ”模型 PICP PINAW ACE CWC 中 的 白 噪 声 剔 出 后 预测 数据 余下 包含 有 效 信息 的 部 分 ， 针 对 
模型 1 0.5116 0.1571 -0.2884 2.8714e+5 剔 出 的 白 噪声 序列 用 Bootstrap 随机 打 乱 后 与 预测 结果 重 构 ， 
模型 2 0.2791 0.1706 -0.5209 3.4983et10 利用 分 位 数 构造 区 间 预 测 结果 。 
ee 模型 3 0.6512 0.1680 -0.1488 286.7808 通过 对 中 国 两 处 机 场 货运 量 ee 发 现 : 四 应 
模型 4 0.7674 0.1694 -0.0326 1.0321 根据 数据 特点 ， 选 择 适合 的 预测 模型 。 通 过 对 航空 货运 需求 
模型 5 0.8605 0.2176 -0.1395 233.3557 的 预测 来 看 , SARIMA 模型 在 季节 分 量 预 测 中 取得 了 比 SVR 
模型 6 0.8372 0.1605 ”0.0372 0.1605 更 好 的 效果 。 针 对 结构 比较 简单 的 时 间 序 列 数据 ，SVR 能 得 
模型 1 0.5581 0.2666 -0.2919 5.8018e+5 到 较 好 的 效果 ,@@ 相 比 于 采用 单一 的 预测 模型 预测 各 个 分 量 ， 
模型 2 0.3256 0.2673 -0.5244 6.5255et10 采用 多 种 预测 模型 分 别 预 测 相 应 的 分 量 ， 充 分 利用 不 同 预测 
六 模型 3 0.6977 0.2680 -0.1523 544.5815 模型 的 天 势 ， 可 以 获得 更 好 的 预测 效果 。SVR-SARIMA 模型 
模型 4 0.8140 0.2488 -0.0360 1.7574 有 效 地 降低 了 单一 模型 对 于 复杂 数据 拟 合 能 力 的 不 足 ， 在 预 
模型 5 0.8605 0.2429 -0.1395 260.4901 测 分 析 中 能 较 好 地 预测 数据 ， 优 于 单一 的 SVR，SARIMA 模 
模型 6 0.8605 0.2669 0.0105 0.2669 型 ， 且 在 预测 方面 比 传统 机 器 学 习 中 的 支持 向 量 机 更 有 效 。 
模型 1 0.6279 0.3234 ”-0.2721 2.6192e+5 名 将 白 噪声 分 量 从 原始 数据 中 剔 出 可 明显 减轻 数据 预测 难度 ， 
模型 2 0.4884 0.3229 -0.4116 2.8015e+8 是 高 预测 精度 。 多 数 原 始 数 据 中 均 包 含有 不 可 预测 的 白 噪声 
0 模型 3 0.8140 0.3221 ” -0.0860 24.1131 即 随机 项 ,将 噪声 分 量 吻 出 即 可 得 到 含有 有 效 信息 的 分 量 ， 
模型 4 0.8837 0.3194 -0.0163 1.0402 此 时 再 进行 预测 便 不 受 随机 项 的 干扰 。 综 上 所 述 ， 所 建立 的 
模型 5 0.8837 0.2783 -0.1163 93.4993 模型 具有 较 好 的 准确 性 、 鲁 棒 性 和 操作 上 的 简单 性 ， 可 以 为 
模型 6 0.9070 0.3238 0.0070 0.3238 区 间 预 测 提供 新 的 思路 。 此 外 ， 鉴 于 所 提出 模型 在 预测 精度 
模型 1 0.7209 0.3835 -0.2291 36140.2084 和 稳定 度 两 方面 的 优势 ， 其 未 来 在 周期 性 数据 预测 等 方向 将 
模型 2 0.8140 0.3857 ” -0.1360 347.4890 有 广阔 的 应 用 前 景 。 
a 模型 3 0.8372 0.3824 -0.1128 107.9575 参考 文献 ; 
模型 4 0.9302 ”0.3771 -0.0198 1.3904 
模型 5 0.9070 0.3306 -0.0930 34.9479 [1] 文 军 ， 藉 由 辉 ， 方 文清 . 航空 货运 量 的 优化 组 合 预测 模型 [四 计算 
模型 6 0.9535 ”0.3842 0.0035 0.3842 机 工程 与 应 用 , 2010, 46 (15): 215-217, 229. (Wen Jun, Jiang Youhui, 
7 区间 预测 结果 评价 (上 海 机 场 ) Fang Wenqing. Optimal combination forecasting model of air cargo 
JI Tab.7 Evaluation of interval prediction results (Shanghai airport) volume [J]. Computer Engineering and Applications, 2010, 46 (15): 215- 
a PINC 模型 ”PICP PINAW ACE CWC 217, 229.) 
模型 1 0.5349 0.1588 -0.2651 90717.9190 [2] 徐 攻 : 薛 锋 . 基于 GM (1, 1) 残 差 二 次 修正 的 铁路 货运 量 预测 [J]. 
模型 2 0.5116 ”0.1578 -0.2884 2.8837e+5 通 运输 工程 与 信息 学 报 , 2019, 17 (02): 44-50. (Xu Li, Xue Feng. 
模型 3 0.6512 0.1585 -0.1488 270.5045 Prediction of rail freight volume of China based on secondary residual 
国人 模型 4 0.7442 0.1534 -0.0558 2.6532 error modification of GM (1, 1) model [J]. Journal of Transportation 
模型 5 0.8837 0.2052 -0.1163 68.9394 Engineering and Information, 2019, 17 (02): 44-50.) 
模型 6 0.8140 0.1580 ”0.0140 0.1580 [3] 瞧 乃 丹 , 向 万 里 , 孟 学 雷 , 等 . 基于 小 波 灰色 GM (1, 1) 模型 的 货运 
模型 1 0.5581 0.1778 -0.2919 3.8700e+5 量 预测 研究 [四 . 铁道 科学 与 工程 学 报 , 2017, 14 (11): 2480-2486. 
模型 2 0.5581 0.1760 -0.2919 3.8288e+5 (Cui Naidan, Xiang Wanli, Meng Xuelei, et al. Railway freight volume 
Be 模型 3 0.6977 0.1768 -0.1523 359.2225 forecasting based on grey GM (1, 1) model and wavelet de-noising [J]. 
模型 4 0.8605 0.1732 0.0105 0.1732 Journal of Railway Science and Engineering, 2017, 14 (11): 2480-2486.) 
模型 5 0.8837 0.2291 -0.1163 76.9557 [4] 游 士兵 ， 严 研 . 逐步 回归 分 析 法 及 其 应 用 [J]. 统计 与 决策 , 2017 
模型 6 0.8605 0.1802 0.0105 0.1802 (14): 31-35. (You Shibing, Yan Yan. Stepwise regression analysis and its 
模型 1 0.6512 0.2316 -0.2488 58641.0122 application [J]. Statistics and Decision, 2017 (14): 31-35.) 
模型 2 0.6744 0.2310 -0.2256 18280.5976 [5] 刘 鞭 ， 赵 静 云 . 基于 系统 动力 学 的 建筑 碳 排放 预测 研究 [J]. 科技 管 
OO 模型 3 0.8372 0.2381 -0.0628 4.7372 理 研究 , 2018, 38 (09): 219-226. (Liu Jing, Zhao Jingyun. Building 
模型 4 0.9070 0.2524 0.0070 0.2524 carbon emissions prediction research based on system dynamics [J]. 
模型 5 0.9302 0.2625 -0.0698 8.8542 Science and Technology Management Research, 2018, 38 (09): 219-226.) 
模型 6 0.9070 0.2372 ”0.0070 0.2372 [6] 刘 笑 你, 任 严 . 基于 不 确定 理论 的 铁路 货运 需求 预测 [J]. 运筹 与 管 
模型 1 0.7209 ”0.3620 -0.2291 34113.3438 理 , 2020, 29 (03): 135-141. (Liu Xiaotong, Ren Shuang. Railway freight 
模型 2 0.8372 0.3746 -0.1128 105.7606 demand forecast based on uncertainty theory [J]. Operations Research 
ee 模型 3 0.8837 0.3783 -0.0663 10.7793 and Management Science, 2020, 29 (03): 135-141.) 
模型 4 0.9302 0.3654 -0.0198 1.3470 [7] 张 亚 男 , 席 洋 , 杨 嘉 钰 ,等 . 海南 自 贸 港 背景 下 海口 市 对 外 客运 量 预 
模型 5 0.9302 0.3118 -0.0698 10.5177 测 研 究 [J]. 交通 运输 系统 工程 与 信息 , 2021, 21 (03): 260-267, 281. 
模型 6 0.9535 0.3899 0.0035 0.3899 (Zhang Yanan, Xi Yang, Yang Jiayu, et al. Forecast of Haikou external 
4 结 束 语 passenger volume considering free trade port in Hainan province [J]. 
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